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ABSTRACT 

The halo region of M31 exhibits a startling level of stellar inhomogeneities, the most prominent of 
which is the "giant southern stream". Our previous analysis indicates that this stream, as well as 
several other observed features, are products of the tidal disruption of a single satellite galaxy with 
stellar mass ~ 10 9 M less than 1 Gyr ago. Here we show that the specific observed morphology of the 
stream and halo debris favors a cold, rotating, disk-like progenitor over a dynamically hot, non-rotating 
one. These observed characteristics include the asymmetric distribution of stars along the stream cross- 
section and its metal-rich core/metal-poor sheath structure. We find that a disk-like progenitor can 
also give rise to arc-like features on the minor axis at certain orbital phases that resemble the recently 
discovered minor-axis "streams", even reproducing the lower observed metallicity of these streams. 
Though interpreted by the discoverers as new, independent tidal streams, our analysis suggests that 
these minor-axis streams may alternatively arise from the progenitor of the giant southern stream. 
Overall, our study points the way to a more complete reconstruction of the stream progenitor and its 
merger with M31, based on the emerging picture that most of the major inhomogeneities observed in 
the M31 halo share a common origin with the giant stream. 

Subject headings: galaxies: M31 - galaxies: interactions - galaxies: kinematics and dynamics 



1. INTRODUCTION 

The relative proximity of the Andromeda galaxy (M31) 
and the global perspective from our external vantage 
point make M31 an excellent laboratory for studying 
the stellar halos of large galaxies. Resolved stellar 
maps of M31's halo, assembled over the past decade, 
have revealed highly complex inhomogeneities, the most 
striking of which is the Giant Southern Stream (GSS), 
extending ~ 150 kpc away from M31 's center in the 
southeast direction ( Ibata et al. 2001; Ferguson et al.l 



120021 : IMcConnachie et all l200l " llbata et al.l (20071 . here- 
after 107) and falling towards M31's cente r with rela- 
tive r a dial velocities as high as ^^Okms -1 (|Ibata et"aTI 
120041 : iGuhathakurta et all 120061 : iKalirai et alT l2QQ6). 
Other significant morphological and kin ematic features 
in th e M31 halo include stellar shelves dFerguson et al.l 
l2002MFardal et all [20071 : iGilbert et all I2QQ7L hereafter 
IF07MG071 ) as w ell as the recently discovered minor axis 
"streams" (|I07l ). The GSS is especially notable because it 
offers an opportunity to precisely measur e M31's poten- 
tial (jlbata et al.ll2004l : iFardal et alJl2QQA hereafter HqI) 
and provides a view into the most significant Local Group 
galaxy disruption in the last Gyr. 

Models detailing the formation of the GSS agree re- 
markably well with most aspects of the observations, 
and suggest the progenitor had a st ellar mass of ~ 
2 x 1O 9 M ([Font et al.ll2006l . Ifq5 lF07h . Our kinematic 
analysis in F07 finds that seemingly unrelated features 
like the "Northeast Shelf" and less prominent "West- 
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ern Shelf" are also the result of the same disruption 
process (|F07f ). a conclusion supp orted by independent 
studies of their stellar p opulations ([Ferguson et al .1 12005: 
iRichardson et al1l2008l ). The observed GSS's most strik- 
ing point of contrast with the models is its asymmetry 
in the transvers e direction. As shown bot h with photo- 
metric samples (IMcConnachie et al.l 12003) and spectro- 
scopic surveys its stellar distribution is sharply 
truncated on the NE side and falls off much more slowly 
on the SW side. In addition, the current models do not 
address t he observed stellar population gradients w ithin 
the GSS ([Ferguson et al.ll2002l : [McConnachiell2006l . 107). 

In this letter, we show that this structure in the GSS 
can be accounted for if the progenitor hosted a cold, ro- 
tating stellar disk, unlike the simple spherical progeni- 
tors used in previous simulations. Surprisingly, we find 
that the disruption of a disk galaxy can also give rise to 
features similar to the recently discovered arc-like minor- 
axis "streams" , leading to the tantalizing possibility that 
most of the major inhomogeneities observed in the M31 
halo are tidal debris from the same galaxy that caused 
the GSS. In Section 2, we briefly describe our model for 
the progenitor and our TV-body study of its tidal dis- 
ruption. In Section 3, we show results from these sim- 
ulations, focusing on the transverse density profile of 
the GSS, the metallicity gradient, and arc-like structures 
that overlap the minor axis. Section 4 summarizes our 
conclusions. 

2. SIMULATION METHOD 

Our simulations are based on the m e thods worke d out 
in o ur earlier papers: iGeehan et all (|2QQ6j ). IFOR and 
F07l . We use the orbit and potential from Table 1 of 



F07I and their spherical Plummer model to represent a 
non-rotating progenitor. For runs with a disk progen- 
itor, we use the same initial position and velocity, but 
substitute a different initial structure of the satellite. 
Briefly, our disk models assume the satellite is com- 
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Fig. 1. — (a): Stellar surface density /met allicity map of M31 from lI07l . The shallower INT/WFC survey is used inside the large ellipse, 
and the deeper CFHT/MegaCam survey outside it. The minor-axis "streams" C and D are visible at the lower left, projecting from the 
larger plume of the GSS. These streams (green) and the GSS cocoon (red) are observed to be more metal poor than the GSS core (yellow), 
(b): Mass surface density map from model Disk A, 840 Myr into the run. The map is 160 kpc on a side, and a dotted contour indicates 
M31's disk orientation. The square indicates the region shown in Figure [2] (c): Same for model Disk B, at 680 Myr. (d): Same for the 
Plummer model, at 840 Myr. (e): Map of the metallicity as a function of position in Disk A (at 840 Myr), with red denoting the highest 
metallicity, dark blue intermediate, and light blue the lowest (see Figure [3ji for a quantitative scale). Boxes indicate the regions used for 
the metallicity histograms in Figure [3] (f): Same for Disk B (at 680 Myr). 



posed of a bulge and rotating disk of stars. For sim- 
plicity we assume that the dark matter associated with 
the galaxy has been tidally stripped before the encounter 
modeled here. We use a hot exponential sech 2 disk with 
mass 1.8 x 10 9 Mq, radial scale length 0.8 kpc, and verti- 
cal scale height 0.4 kpc. We add to this a Hernquist bulge 
of mass 4 x 10 8 M and scale length 0.4 kpc. We initial- 
ize both components with the package ZENO written by 
Josh Barnes. We evolve the satellite in M31's potential 
starting from 12 evenly spaced orientations of the disk. 
From this we select two models displaying particularly 
good agreement with observational features, referred to 
here as Disks A and B. In a forthcoming paper we will 
conduct a more systematic survey of possible initial con- 
ditions and quantify the debris structure in detail (Fardal 
et al., in preparation). 

The remaining det ails o f the simulations are the same 
as those in lF06l and lF07l . We set the satellite in mo- 
tion inbound and slightly past apocenter, minimizing ini- 
tial transients from M31's tidal forces. We run the sim- 
ulations with the multistepping tree code PKDGRAV 
([Stadell 12001). Our simulations include the self-gravity 
of the progenitor satellite, but ignore dynamical friction, 
perturbations from the other M31 satellites, and the his- 
tory of the progenitor prior to the orbit that produces 
the GSS, as justified in our earlier papers. We stop the 



simulations 840 Myr into the run, which is approximately 
65 My r past the initial pericenter. This time was chosen 
in lF07l to give a reasonable match to the debris structure 
around M31, including the radii of the "shells" on the E 
and W sides. 

3. RESULTS 
3.1. Stream morphology 

Figures [Tj3,d show surface density maps based on the 
Disk A and Plummer models, respectively. Both models 
reproduce the main feature of a stream extending to the 
SE. They also reproduce the observed line-of-sight dis- 
tances and velocities along the GSS. However, the trans- 
verse distribution of GSS stars is strikingly different be- 
tween the two models — Disk A displays a much sharper 
NE edge. The observed star-count maps ([Ferguson et al.l 
2002, lIQTh are not directly comparable since they con- 
tain both non-GSS-related M31 components and non- 
M31 contaminants and are not explicitly calibrated to 
stellar surface density, but the morphology of the GSS in 
these maps appears much closer to our disk model. 

Figure [T]i shows that the Plummer model results in 
a large amount of stars spilling over as far as th e SE 
minor axis, located to the NE of the GSS. When [GOT! 
compared their Keck/DEIMOS spectroscopic data near 
M31's SE minor axis to this model, they noted much 
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Fig. 2. — Comparison of the minor-axis contamination to the 
observations of G07. The IG07I DEIMOS masks (rectangles) are 
grouped into inner minor-axis masks, outer minor-axis masks, and 
a single mask (fl35) offset from the minor axis. M31's center is 
at (0, 0). The inset plots for each group show the ratio Rm of the 
strength of the GSS component to the peak of the GSS at the same 
-Rproj • Rm is measured as discussed in the text for Disks A and B 
and the Plummer model 840 Myr into the runs. The observational 
estimates and ±2cr error bars from Gilbert et al. (in preparation) 
are plotted as horizontal solid and dotted lines. The Plummer 
model clearly contributes too much debris on the minor axis. 



less spillover from the GSS than predicted by the model. 
Gilbert et al. (in preparation) has quantified this by di- 
viding the number of stars moving with GSS-like veloc- 
ities on the minor axis to those in the GSS core at the 



same projected radius R pr0 y For the nine innermost 
DEIMOS masks on the minor axis combined, this ra- 
tio Rm = 0.01 ± 0.02; for the three outermost masks 
on the minor axis, R m = 0.05 ± 0.02; and for the mask 
fl35 located somewhat nearer the GSS, they find a likely 
detection of GSS material with Rm = 0.17 ± 0.06. 

In Figure [2j we compare the density of GSS stars in all 
three TV-body models to these results. We have selected 
"GSS" particles by defining the trend of radial veloc- 
ity v with i? pro j and then taking stars that fall within 
±80kms _1 of this velocity in the given field. We also re- 
strict the particles to those actually in the GSS's "shell". 
We then repeat the procedure for a control field located 
at the peak of the GSS at the same i? pro j, using a smaller 
interval ±40kms _1 as the GSS core has a sharper veloc- 
ity distribution. Clearly the two disk models are in better 
accord with the observations than the Plummer model. 

The sharper NE edge and smaller minor-axis contami- 
nation of the disk models thus imply that the progenitor 
was rapidly rotating. We will explore this argument in 
more detail in Fardal et al. (2008, in preparation). 

3.2. Metallicity pattern 

The mean color of GSS RGB stars is observed to 
vary in the transverse direction: t he GSS is significantly 
broader in blue tha n in red stars ([Ferguson et al.l [20021 : 
lMcConnachiell2006t) . This is probably due to a metallic- 
ity gradient. 1071 quantified the metallicity distribution 
in two GSS-dominated regions, one in the center of the 



GSS and one in a less dense "cocoon" region to the SW, 
and showed that the latter has a lower mean metallicity. 

Disk galaxies, of course, tend to have metallicity gra- 
dients. Therefore it is interesting to see how a plausible 
gradient in our disk progenitor translates to the metal- 
licity pattern on the sky. 

We use a simple parametric model to produce a plau- 
sible metallicity gradient in our initial disk model. We 
first find the specific orbital energy Ei of each parti- 
cle. We then assign it a metallicity using [Fe/H] = 
A z + B z log 10 [-^/(SO kms -1 ) 2 ], setting A z and B z 
to agree with the results of Il07l as explained below. This 
produces the metallicity gradient seen in Figure [3^. Ob- 
servational results for the stars in the small disk galaxy 
M33 are also plotted, with the radius for both galaxies 
normalized by the disk scale length; the metallicity pat- 
tern of our disk model agrees quite well. The GSS pro- 
genitor should perhaps be lower in metallicity than M33 
by a few tenths of dex due to its lower inferred mass, 
but the photometric metallicity measurements probably 
have systematic uncertainties at this level in any case. 

Figure []J shows the sky view of the resulting model 
metallicity pattern. The gradient along the stream is 
very weak, but the mean metallicity along the denser 
central part is clearly higher than in the broad wing 
to the SW , similar to the pattern seen in M31's GSS. 
Using Il07l 's Figure 27, we estimate the "core" and "co- 
coon" regions (at i? pro j ~ 60 kpc) have mean metallic- 
ities of [Fe/H] = —0.54 and —0.71, respectively (mean 
[Fe/H] = -0.51 was obtained at the GSS' sharp NE 
edge by iGuhathakurta et "all l20Q6h . Figure []J shows 
"broad wing" and "central GSS" boxes chosen at a sim- 
ilar radius, but better matching the slightly different 
model stream position. Once we set Az = —0.70 and 
Bz = 1.06, the metallicities in these boxes are also —0.54 
and —0.71. The bare fact we can match two metallici- 
ties with two parameters is not in itself meaningful, but 
it is significant that the magnitude and sign of our ini- 
tial metal gradient are very reasonable (Fig. [3^). Fig- 
ure Od shows that within each box there is a wide range 
of metallicities; the distributions in lI07l appear somewhat 
broader, but given measurement errors and the contribu- 
tions from other halo components this is not surprising. 

3.3. Minor- axis arcs 

Using their MegaCam photometric survey of M31's 
halo, 107 found multiple surface density ridges along the 
minor axis which they called "streams" . Streams C and 
D (the two closest to M31) form a pair of curving ridges 
at slightly different orientations, which appear to merge 
as they approach the survey boundary (see their Fig. 22). 
Stream C appears to be slightly broader than stream D, 
and slightly more metal-rich, though not as metal-rich as 
the GSS core/cocoon. From I07's Figure 33 we estimate 
the mean metallicity of streams C and D to be —0.82 
and —0.91 respectively. iMori fc Rich! (|2008l ) suggested 
these "streams" might be shell features from a satellite 
disruption, similar to the event that created the GSS but 
from a different progenitor. 

While studying our overall sample of runs based on 
12 disk orientations, we noticed one (Disk B) contain- 
ing two curious "arcs" crossing the minor axis. These 
arcs are clearest at the step 680 Myr into the run shown 
in Figure [Tfc. Morphologically, the two arcs somewhat 
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Fig. 3. — (a): Metallicity values in the original disk prior to disruption are shown by black dots, wher e the radius is plotted in u nits of disk 
scale length. For comparison, obs erved results for M 33's disk stars are plotted as symb ols and lines: [Stephens & Frogel (2002) (square); 
linear approx i matio n to points of Ki m et aTl (|2002T ) (straight line); Galleti et al. (2004) (diamonds); McConnachie et al. (2006) (cross); 
Barker et al. (2007) (triangles). The colorbar translates [Fe/H] to the color scale of Figures [T^f. (b): Histogram of particle metallicity 
values in model Disk A within the "core" and "cocoon" regions marked in Figure ^.(c): Histogram of particle metallicity values in model 
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Disk B within the "core", "cocoon",. "N arc" and "S arc" regions marked in Figure [If. . 

resemble streams C and D, with a tatter southern arc In summary, several strands oi observational evidence 



nearly merging into a sharper northern arc. Like the ob- 
served "streams" , neither arc crosses the GSS to the SW. 
Compared to the observed arcs, the simulated arcs are 
significantly further from M31's center. 

As Figure [TJ shows, the simulated arcs are significantly 
less rich in metals than the GSS. Using the same metal- 
licity model as for Disk A and the regions defined by 
boxes in this figure, the mean [Fe/H] is —0.78 for the 
southern arc and —0.90 for the northern arc. Thus there 
is considerable if inconclusive evidence that these arcs 
are close analogues of the "streams" in Il07l . 

In our model, these two arcs originate from the outer 
regions of the disk, and are sharp mainly because of the 
relatively cold velocity field of the disk. Both arcs consist 
of material that takes a path around M31's center nearly 
opposite to the bulk of the progenitor, explaining why 
they lie so far from the GSS. The large size of our disk 
is thus crucial; a compact progenitor resembling M32, 
for example, would be unable to produce similar arcs. 
The southern arc consists of a group of particles sharing 
nearly the same energy, and come from fairly far out 
in the progenitor's disk. The northern arc consists of 
particles that lie even further out (explaining its lower 
metallicity on average), which form a tidal tail during 
the interaction with M31. 

We cannot yet explore the full parameter space of the 
collision for the presence and properties of these arc-like 
features. However, we did conduct a few additional runs 
with changes to the disk mass, radius, and orientation 
of Disk B, finding the arcs were sensitive to the exact 
input parameters. Thus we will require more theoretical 
investigation as well as more observational constraints 
to determine whether the arcs explain some of the Il07l 
minor- axis streams, or are merely a fortuitous similar- 
ity. If the arcs are shown to be related to the GSS, they 
will be a very solid argument for the disk nature of the 
progenitor, and will place strong constraints on the pa- 
rameters of the collision. 

4. CONCLUSIONS 



suggest that the GSS originated from a progenitor with 
a strong sense of rotation, such as a disk galaxy. The 
transverse density profile of the GSS is more easily pro- 
duced by a rotating satellite. The observed decline in 
mean metallicity from the central core of the GSS to its 
"cocoon" to the SW suggests that the progenitor had 
a strong radial metallicity gradient, of the sort found 
mainly in disk galaxies. Furthermore, several observed 
arcs lying across the minor axis in M31 have very sugges- 
tive analogues in one of our runs. If shown to be related 
to the GSS in the manner suggested by our model, these 
features would be definite confirmation of a disk-like pro- 
genitor. 

The notion of a disk galaxy progenitor is somewhat at 
odds with age measurements of the GSS, w hich suggests 
little s tar formation during the last 4 Gyr (| Brown et al.l 
2006a.b). However, the fields used to infer this were 
placed in the central, metal-rich part of the GSS; it is 
possible that the progenitor had an age gradient as well 
as a metallicity gradient, with the older stars on the in- 
side. Age measurements in the GSS cocoon would there- 
fore be interesting. It is also possible that the GSS pro- 
genitor was more similar to an SO galaxy than a spiral, 
perhaps due to stripping of its gas in an earlier phase of 
its encounter with M31. 

Many papers have used metallicity to assess the rela- 
tionship among various M31 disk and halo features. Our 
suggestion that the GSS progenitor had a strong metal- 
licity gradient means that metallicity can no longer be 
used as a reliable fingerprint of origin. This obviously 
complicates the forensic reconstruction of M31's merger 
history. Despite this, the rapidly growing database 
on M31 halo structure is a fascinating puzzle, offering 
unique insight into the life of a typical disk galaxy and 
the death of its unfortunate former companions. 

We thank Tom Quinn and Joachim Stadel for the use 
of PKDGRAV, Josh Barnes for the use of ZENO, and 
Alan McConnachie and Roger Davies for helpful conver- 
sations. 
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